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Abstract: The electromagnetic and gravitational quasinormal spectra of (3 + 1)- 
dimensional plane-symmetric anti-de Sitter black holes are analyzed in the context of the 
AdS/CFT correspondence. According to such a correspondence, the electromagnetic and 
gravitational quasinormal frequencies of these black holes are associated respectively to the 
poles of retarded correlation functions of /^-symmetry currents and stress-energy tensor in 
the holographically dual conformal field theory: the (2 + l)-dimensional = 8 super- Yang- 
Mills theory. The connection between AdS black holes and the corresponding field theory 
is used to unambiguously fix the boundary conditions that enter the proper definition of 
quasinormal modes. Such a procedure also helps one to decide, among the various different 
possibilities, what are the appropriate gauge- invariant quantities one should use in order 
to correctly describe the electromagnetic and gravitational blackhole perturbations. These 
choices imply in different dispersion relations for the quasinormal modes when compared 
to some of the results in the literature. In particular, the long-distance, low-frequency limit 
of dispersion relations presents the characteristic hydrodynamic behavior of a conformal 
field theory with the presence of diffusion, shear, and sound wave modes. There is also 
a family of purely damped electromagnetic modes which tend to the bosonic Matsubara 
frequencies in the long- wavelength regime. 
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1. Introduction 



1.1 Motivations and overview 

Theoretical studies on black holes in asymptotically anti-de Sitter spacetimes have 
attracted substantial attention since the advent of the anti-de Sitter /conformal field theory 
(AdS/CFT) correspondence [1-3]. In particular, the quasinormal-mode (QNM) spectra 
of various types of asymptotically AdS black holes have been analyzed since then (see 
Refs. [4-30] for a sample). According to the AdS/CFT correspondence, an asymptotically 
AdS black hole is, in the CFT side, associated to a system in thermal equilibrium whose 
temperature is the Hawking temperature of the black hole. In such a context, blackhole 
perturbations correspond to small deviations from equilibrium of the CFT thermal system, 
and the characteristic damping time of perturbations, which is given by the inverse of the 
imaginary part of the fundamental QNM frequency, is a measure of the dynamical timescale 
of approach to thermal equilibrium of the corresponding conformal field theory [31]. 

The literature on QNM of AdS black holes includes studies taking into account a vari- 
ety of different aspects such as the topology of the event horizon, the number of dimensions 
of the spacetime, the particular type of perturbation fields considered, and also the special 
parameters which characterize each different black hole itself. Each one of these variant 
properties rcficcts on the dual CFT. For instance, assuming the (3 + l)-dinicnsional AdS 
spacetime contains a plane-symmetric black hole, then the holographic field theory is de- 
fined over the (2 -|- l)-dimensional Minkowski spacetime, which is the conformal boundary 
of the bulk AdS spacetime. Moreover, different blackhole parameters characterize differ- 
ent dual plasmas in the CFT side, and different equilibrium states of such systems at the 
boundary. 

An important issue in the study of the vibrational modes of black holes is the choice of 
appropriate boundary conditions. In the case of asymptotically flat spacetimes, the solu- 
tions to the wave equations governing linear perturbations are, near the boundaries, given 
by plane wave functions. QNM are then defined as solutions which satisfy physically well 
motivated boundary conditions, namely, purely ingoing waves at the horizon and purely 
outgoing waves at infinity (see Refs. [32,33] for reviews). For anti-de Sitter black holes, on 
the other hand, the condition at the future horizon is the same as for asymptotically flat 
spacetimes, but now there are no natural conditions to be imposed on the perturbation vari- 
ables at the AdS infinity. These can be Dirichlet, Neumann, or Robin boundary conditions, 
depending on whether it is required that the field perturbations, their derivatives or a com- 
bination of both vanish at the AdS boundary, respectively. In the study of the evolution of 
a massless scalar field in (3-M)-, (4-1-1)-, and (6-|-l)-dimensional Schwarzschild-AdS space- 
times, Horowitz and Hubeny [31] computed the corresponding quasinormal-mode spectra 
by imposing Dirichlet boundary conditions on such a field at infinity. This option was well 
justified in that context, since by writing the radial part of the Klein- Gordon equation in a 
Schrodinger-like form, the resulting effective potential diverges at that boundary. The same 
boundary condition was used to study massless scalar and electromagnetic perturbations 
of (2 -|- l)-dimcnsional Baiiados-Teitelboim-Zanelli (BTZ) black holes [34]. For BTZ black 
holes, an analytical closed form for the quasinormal frequencies was derived [35], and it 
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was verified that the quasinormal frequencies correspond exactly to the poles of retarded 
correlation functions in the dual (1 + l)-dimensional CFT [36]. It was also suggested in 
Ref. [37] that the relation between quasinormal modes and singularities of correlation func- 
tions should also hold for scalar fields in higher-dimensions, as far as the frequencies are 
computed by imposing Dirichlet boundary conditions on such fields at AdS infinity. 

In the meantime, two fundamental difficulties arise when considering gravitational 
and/or electromagnetic perturbations of AdS black holes, particularly in higher dimensional 
spacetimes. The first problem is related to the arbitrariness in the choice of gauge- invariant 
perturbation fields. In fact, there is an infinity of gauge-invariant combinations of metric 
(or vector potential) fluctuations that can be used as fundamental variables governing 
the gravitational (or electromagnetic) perturbations. The second problem is related to 
the ambiguity in defining appropriate boundary conditions for the quasinormal modes. 
A traditional way to face such arbitrariness is opting for master variables that lead to 
equations generalizing those for perturbations in asymptotically flat spacetimes. That is to 
say, variables are chosen in such a way to put the radial part of the fundamental equations 
into a Schrodinger-like form. Prom now on, the corresponding master variables shall be 
called the Rcggc-Wheeler-Zerilli (RWZ) variables.^ With such a choice of variables, it 
was investigated gravitational and/or electromagnetic perturbations of the Schwarzschild- 
AdS [38,41-46], Rcissncr-Nordstrom-AdS [47-49], and Kerr- AdS [50] black holes, as well as 
the perturbations of black holes with non-spherical topologies [51,52], including the plane- 
symmetric ones [53-55] . Analogously to the massless scalar field case, in all of these works 
the quasinormal modes were computed by imposing Dirichlet boundary conditions on the 
master fields at infinity. Alternative boundary conditions for the same Regge-Wheeler- 
Zerilli variables have been discussed in Refs. [56,57]. 

A different route was taken by Nunez and Starinets [58] , who defined the quasinormal 
frequencies of a perturbation in an asymptotically AdS spacetime as "the locations in the 
complex frequency plane of the poles of the retarded correlator of the operators dual to 
that perturbation". To compute the real-time correlation functions, they suggested using 
the Lorentzian AdS/CFT prescription of Refs. [37,59]. The quasinormal-mode definition 
supplied by Nunez and Starinets was explored in Ref. [60] , where a new set of fundamental 
variables was introduced to study electromagnetic and gravitational perturbations of (4 -|- 
l)-dimensional plane-symmetric black holes (or black branes, for short). It was shown 
there that the imposition of Dirichlet boundary conditions on such a new set of gauge- 
invariant variables at infinity leads exactly to the quasinormal frequencies associated to 
the corresponding black branes. In the present work these kind of fundamental variables 
shall be called the Kovtun-Starinets (KS) variables. 

An important consequence of the Nuhez-Starinets approach [58] is that the resulting 
quasinormal-mode spectra present a set of dispersion relations, here called hydrodynamic 
QNM, that behave like diffusion, shear, and sound wave modes in the long-wavelength, low- 
frequency limit [60]. These results are totally consistent with what is expected from the 

^In the first study of gravitational QNM in AdS spacetimes, Cardoso and Lemos [38] used the same 
kind of variables as the early works in asymptotically flat spacetimes by Regge and Wheeler [39], and by 
Zerilli [40]. 
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CFT point of view, and they provide a non-trivial test of the AdS/CFT correspondence. It 
is also worth noticing that neither the electromagnetic diffusion mode nor the gravitational 
sound wave mode are obtained by imposing Dirichlet boundary conditions on the RWZ 
master variables. For Schwarzschild-AdS and topological- AdS (3 + l)-dimensionaI black 
holes, it was only possible to obtain sound wave modes in the gravitational quasinormal 
spectra by requiring that a specific combination of the master field and its derivative 
vanishes at infinity [57,61,62]. 

1.2 The present work 

1.2.1 General procedure 

In this work the definition of QNM given by Niihez and Starinets [58] is applied to 
compute the quasinormal frequencies associated to electromagnetic and gravitational per- 
turbations of (3 -I- l)-dimensional plane-symmetric AdS black holes. The overall procedure 
is similar to that of Ref. [60] and consists of the following steps: 

(1) Initially the translation invariance of the static plane-symmetric AdS spacetimes is 
used to Fourier transform the fluctuation fields with respect to time and to the two 
Cartesian coordinates {x,y) of the plane. 

(2) With the spatial wave vector chosen to be in the y-direction, both the electromagnetic 
and the gravitational perturbation fields are separated into two sets according to 
their behavior under the transformation x — ^ —x: odd (axial, or transverse), and 
even (polar, or longitudinal) perturbations. 

(3) Each sector of perturbation fields is governed by a set of linearized differential equa- 
tions. In all of the cases studied here, the complete set of perturbation equations can 

be decoupled in order to obtain a unique second-order differential equation, which is 
the fundamental equation of that perturbation sector. The fundamental equations are 
written in terms of gauge-invariant combinations of the perturbation fields, extend- 
ing the original definitions of Kovtun-Starinets variables [60] to (3 -|- l)-dimensional 
spacetimes. 

(4) Then, the standard AdS/CFT prescription of Ref. [37] is applied to express the real- 
time ii-symmetry current and stress-energy tensor correlators in terms of quantities 
which represent the asymptotic behavior of perturbations near the AdS-space bound- 
ary. Such a procedure shows that the imposition of Dirichlet boundary conditions 
on Kovtun-Starinets variables at infinity leads to the poles of the CFT correlation 
functions, and therefore, according to the Niihez-Starinets definition of QNM, to the 
quasinormal spectra of the plane-symmetric AdS black holes. 

(5) With well defined boundary conditions and a set of decoupled fundamental equations, 
the hydrodynamical limit of the QNM spectra is then analyzed. This limit is reached 
for perturbation modes in which the frequency and the wavenumber are much smaller 
than the Hawking temperature of the black hole. 
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(6) The last step is numerically compute the electromagnetic and gravitational quasinor- 
mal dispersion relations for different blackhole parameters. For such a purpose, the 
Horowitz-Hubeny method [31], which reduces the problem of finding QNM frequen- 
cies to that of obtaining the roots of infinite polynomial equations, is used. 

1.2.2 Main results 

Among the new results found in the present work, it is worth mentioning the following 
ones. 

• First, the derivation of the electromagnetic diffusion mode and the gravitational 
sound wave mode is performed by means of a traditional QNM calculation. These 
modes were earlier obtained by Herzog [63,64], who utilized the AdS/CFT prescrip- 
tion [37] to directly compute the hydrodynamic limit of the CFT i?-symmetry current 
and stress-energy tensor correlators. 

• Second, it is found that the procedure of imposing Dirichlet boundary conditions 
on the gauge-invariant KS variables breaks the isospectrality between the axial and 
polar electromagnetic QNM, that follows from RWZ variables. As a result, the po- 
lar electromagnetic quasinormal modes are totally new, since the KS variable with 
Dirichlet boundary condition yields a different spectrum when compared to the RWZ 
variable with the same kind of boundary conditions. Regarding to the electromag- 
netic axial perturbations, the dispersion relations found here enlarge previous results 
of Ref. [53]. 

• Third, it is found in addition that the complete spectra of electromagnetic QNM 
present a tower of purely damped modes which tend to Matsubara frequencies char- 
acteristic to bosonic systems in the long-wavelength regime. However, these quasi- 
normal modes do not exist for all wavenumbers. In fact, there is a saturation value for 
the wavenumber above which the electromagnetic purely damped modes disappear. 

• Fourth, another result to be mentioned is the difference between the spectrum of 
the gravitational polar perturbations, computed by using the KS variable, and that 
obtained using the RWZ master variable [54]. The differences are specially signifi- 
cant when the fluctuation wavenumber is of the same order of the magnitude of the 
blackhole temperature. 

• And last but not least, the dispersion relations calculated here complete the previous 
results for axial gravitational QNM of (3 + l)-dimensional plane-symmetric AdS black 
holes [53,54]. 

1.2.3 Structure of the paper 

The layout of the present article is as follows. Sect. § contains a brief summary of the 
relation between the plane-symmetric AdS4 black holes and the eleven-dimensional super- 
gravity solution associated with a stack of M2-branes. In the sequence a detailed study 
of the electromagnetic quasinormal modes is performed (Sect. ^): The basic equations 
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are obtained in Sect. 3J- and the connection between the blackhole perturbations and the 
CFT i?-symmetry currents is explored in Sect. 3.2; the hydrodynamic modes of the elec- 
tromagnetic perturbations are studied in Sect. 3.4, and the general dispersion relations of 
electromagnetic QNM are reported in Sect. 3.5. The gravitational quasinormal modes are 
studied in Sect. The basic equations are obtained in Sect. 4J; O is devoted to investi- 
gate the relation between the gravitational QNM and the stress-energy tensor correlators 
in the holographic CFT; the hydrodynamic modes of the gravitational perturbations are 
studied in Sect. iA, and the numerical results for the dispersion relations of the remaining 
gravitational QNM are presented in Sect. 4^. The article is completed, in Sect. with 
the analysis and interpretation of the main results. 



1.3 Notation and conventions 

Natural units are going to be used throughout this paper, i.e., the speed of light 
c, Boltzmann constant ks, and Planck constant h are all set to unity, c = ks = ^ = 1- 
Regarding to notation, capital Latin indices M, N, ... vary over the coordinates of the whole 
AdS spacetime, while Greek indices fi, ... label different coordinates at the boundary, 
and small Latin indices i, j, ... vary only over the spacelike coordinates at the boundary. 
The convention for the metric signature and for all the definitions of curvature tensors 
follow Ref. [65]. 



2. M2-branes and the plane-symmetric black holes 

2.1 The background spacetime 

Since the QNM definition of Nuhez-Starinets [58], that is adopted in this work, makes 
heavy use of the relation between AdS black holes and conformal field theories at finite 
temperature, it becomes important to review here how the plane-symmetric AdS4 black 
holes arise in the context of the AdS/CFT conjecture.^ 

A fundamental role in the AdS/CFT correspondence^ is played by extended two- 
dimensional objects known as M2-branes [68]. The world- volume theory of M2-branes 
is a (2 + l)-dimensional non-Abelian Yang-Mills theory which presents = 8 supersym- 
metries in addition to a SU (N) gauge group. The coupling constant of the theory flows 
to strong coupling in the infrared limit, and it is believed that the flow is to an infrared- 
stable fixed point that describes a superconformal field theory [69]. This CFT also has an 
emerging i?-charge symmetry which is expanded to SO (8). 

From the supergravity point of view, a stack of N M2-branes is described by a 
nonextremal solution to the supergravity equations of motion, characterized by the met- 
ric [63, 70, 71] 

^The brief summary presented in this section is based on material found in Refs. [63,66,67]. 
^The interest here is the AdS4/CFT3 correspondence. 
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Hi?) = 1+1^] and = 1 - ) , (2.2) 



where 

Hi?) = l+( 

and by a four-form field whose dual Hodge is given by 

= F7 = 6i2^Vol(5^)e, (2.3) 

where e stands for the Levi-Civita tensor on . According to the AdS/CFT correspon- 
dence [1-3], the (2 + l)-dimensional AA = 8 CFT is dual to M-theory on the background 
spacetime (p.l|). Furthermore, the quantization condition on the -F4 flux connects the 
parameter i? to the number of branes N [72]: 

i?V = A^=^/2^?i^/2, (2.4) 

where k\\ is the gravitational coupling strength in (10 + l)-dimensional supergravity. 

In the large limit (A 3> 1), one can consider only the near-horizon region (r <C K) 
of the spacetime (|2.1|). Function i^(r) then reduces to ff(F) = /r^ . Moreover, defining 



a new radial coordinate by r = r^/2i?, metric ( |2.1| ) becomes 

ds^ = ^ [-W)dt' + dx^ + dy^] + 5^ + R'd^l (2.5) 



The AdS part of the metric (p^), associated to the coordinates {t,j;,y,r}, is identical 
to the solution of Einstein equations with negative cosmological term corresponding to a 
(3 + l)-dimensional plane-symmetric AdS black hole [73-76]: 

2 

ds^ = -fir)dt^ + fir)-^dr^ + ^(dx^ + dy^), (2.6) 
where the horizon function /(r) is given by 

and the seven-sphere radius R has been rewritten as ii = 2L, with L now representing 
the AdS radius of the spacetime (|2.6D . Parameters ro and L are related to the blackhole 
Hawking temperature T by 



2.2 Normalization of the field action 

The full theory is the eleven-dimensional supergravity on AdS4 x S^, and the existence 
of a compact seven-sphere enables one to consistently reduce the theory to Einstein-Maxwell 
theory on AdS4 [67,77,78]. The main objective in summarizing such a procedure here is 
to make explicit the dependence of the action for the fields in the AdS4 spacetime on the 
number of colours N, which is one of the parameters characterizing the holographic CFT. 
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Upon Kaluza-Klein dimensional reduction, the Maxwell gauge field A^j arises from a 
combination of metric and F4 form perturbations in the eleven-dimensional supergravity. 
This field corresponds to a U{1) subgroup of the 50(8) symmetry group of the com- 
plete spacetime ( |2.1| ). The mechanism of dimensional reduction also furnishes the (3 -|- 1)- 
dimensional Einstein-Maxwell action with a negative cosmological constant A = — 

S=^J d'^V^ + ^'^'^'^^""') ' (2-9) 

where TZ denotes the Ricci scalar and Fmn is the electromagnetic strength tensor, and for 
the purposes of the present analysis the electromagnetic Lag rangian ^Cem ^ -^mn-^ is 
considered as a perturbation on the gravitational Lagrangian Cgr ~ TZ+Q/LP'. It is assumed 
that the gravitational coupling constants in four and eleven dimensions are related by means 
of the seven-sphere volume [67], 



1 R'\o\{S'' 



(2.10) 



Then, considering that the volume of a unitary seven-sphere is Vol(5'^) = vr^/3, and using 
the standard normalization ( |2.4D for kh, it is found 

J--^?^ (211) 

Action ( p. 91 ) with the gravitational constant K4 given in terms of the number of colours N 
and of the anti-de Sitter radius L is the desired result, which is needed for the development 
of the present work. 



3. Electromagnetic quasinormal modes 

3.1 Perturbation equations 

In the AdS/CFT context, the electromagnetic field in the AdS bulk couples to the CFT 
i?-symmetry currents at the spacetime boundary. Hence, in order to construct the current- 
current two-point correlation functions in the CFT, it is necessary to consider fluctuations 
of the gauge field A^- Such a field is implicitly defined by 

Fmn — dufAf^ — 9jv^M) (3-1) 



with Fmn satisfying equations of motion derived from the action ( |2.9|) . Therefore, consid- 
ering the electromagnetic field as a perturbation on the background spacetime of metric 
( |2.6| ) , the resulting equations of motion for Am are the usual Maxwell equations 

Om {V^g"^9''^F_,,) = 0, (3.2) 



where g^^ stands for the metric components given by (2.6). 
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When looking for solutions to Eqs. (p.2|), by taking into account the isometrics of the 
background metric ( |2.6| ) , it is convenient to decompose the gauge field in terms of Fourier 
transforms as follows 

A,,{t,^,y,r) = j du:dk^dkye-'^'+'^-^+"'yyA,,{oo,k^,ky,r). (3.3) 

Furthermore, without loss of generality, in the plane-symmetric background spacetime ( ^.(^ ) 
one may choose the wave three-vector k in the form k^ = {ko,kx,ky) = {—Lo,0,q). This 
is carried out through an appropriate rotation in the x — y plane, in such a way that the 
Fourier modes of the gauge field propagate along the y direction only. With such a choice, 
the electromagnetic perturbations Am can be split into two independent sets according to 
their behavior under parity operation, x — > —x: 

• Axial (odd, or transverse) perturbations: Ax', 

• Polar (even, or longitudinal) perturbations: At, Ay, A^. 

Since these two sets of perturbations are orthogonal sets, they can be studied separately, 
as it is done in the following. 

3.1.1 Equations for axial perturbations 

Axial electromagnetic perturbations are governed by the transverse component of 
Maxwell equations ( p.2D , which gives 

^dPAx ^ dfdAx ^ / c^V^-g^LV \ ^ _^ 
dr"^ dr dr \ fr"^ J ^ ' 

where, to simplify notation, the tilde was dropped. Ax — > Ax- Moreover, it follows from 
Eqs. ( PI ) and (^^), together with k^ = (— w, 0, g), that Ax is proportional to the transverse 
component of the electric field: Ex = iujAx- Therefore, being a gauge-invariant quantity. 
Ax is also a good candidate as master variable for axial perturbations. In fact, it is possible 
to cast Eq. (|3.4| ) into a Schrodinger-like form [53] 



dr^ I \ r 

where ^*~^(r) = Ax{r), and the tortoise coordinate r* is defined in terms of the radial 
coordinate r by 

^ = fir)- (3.6) 
dr^ 

For the present purposes it is convenient to change coordinates to the inverse radius 
u = ro/r, and then, by writing Eq. (l3.5|) in terms of Ex it results 

< + |< + "^^Ex = 0, (3.7) 
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where the primes denote derivation with respect to the variable u, and tr and q are respec- 
tively the normalized frequency and wavenumber, defined by 



3a; , 3(? , , 

"' = 4^ q = ^, (3.8) 



with T being the Hawking temperature of the black hole, given by Eq. (p.^). Function f) 
is obtained from Eqs. ( |2.2[) , or from Eq. ( |2.7] ), and in terms of the variable u = tq/t reads 



i) = \){u) = l-u\ (3.9) 

Quantity Ex shall be the fundamental gauge-invariant variable to be used in the present 
analysis of the QNM modes for axial electromagnetic perturbations. 

3.1.2 Equations for polar perturbations 

Differently from the axial electromagnetic perturbations, the components of the gauge 
field Am corresponding to the set of polar perturbations are not gauge invariant. The gauge 
freedom can then be used in order to simplify the relevant equations of motion. In fact, the 
invariance of Maxwell equations under the gauge transformation A^ + allows 

choosing A in such a way that one of the components At, Ay, or Ar vanishes. For instance, 
it is possible to work in the so-called radial gauge, in which Ar = [63]. In this gauge. 
Maxwell equations ( |3.2[) corresponding to the polar perturbations are 

ujr^±At + qL'f^Ay = 0, (3.10) 
ar ar 

^'-^At + 2r^At - J {qu^Ay + q^A^) = 0, (3.11) 

Notice that this set of equations does not constitute a linearly independent system, and 
any subset composed by two of such equations determines the two remaining unknown 
components of the gauge field. At and Ay. 

From now on, At or Ay could be adopted as a primary variable and equations ( ^.10 )- 



( 3.12 ) could be decoupled in order to find a unique differential equation for one of these 
functions. However, to avoid the inconvenient of dealing with gauge dependent quantities, 
it is interesting to use the electric field components, which are gauge- invariant quantities. 
Even though this choice eliminates gauge-dependent potential fields, a residual ambiguity 
is left: from the electric field components Er = dAt/dr and Ey = i{qAt + ooAy), what is 
the best choice? 

The answer to the last question is not simple and both of the possible answers have 
been tried in the literature. For instance, inspired by preceding works studying blackhole 
perturbations in asymptotically fiat spacetimes [79], Cardoso and Lemos opted for the 
radial component E^ in studying QNM of Schwarzschild-AdS black holes [38], and plane- 
symmetric AdS black holes [53]. Introducing a new variable ^'*+'(r) = r'^Er{r) and using 
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Eqs. ( ^l0|) and they were able to reduce the system of equations into a unique 

ordinary differential equation of Schrodinger type 



(|^ + .^)^(-) = /(^)^«, (3.13) 

which has the same form as the fundamental equation for axial perturbations, Eq. (p.5|). 
An interesting consequence of this fact is that the QNM spectra for both the axial and 
polar perturbations, with the same boundary conditions, are identical [38,53]. As a matter 
of fact, an open question left behind in the early works computing electromagnetic QNM of 
AdS black holes is the lack of a convincing physical justification for the choice of Dirichlet 
boundary conditions at infinity for both the polar and the axial perturbations. Such an 
issue will be considered in the sequence of this work. 

As far as one is interested in computing the polar electromagnetic quasinormal frequen- 
cies, it will be shown in Sect. that Ey is more appropriate as a fundamental variable 
than Er- This was the choice made, for instance, by Kovtun and Starinets in studying 
QNM of black branes in (4 + l)-dimensional spacetimes [60]. Following these authors, in 
the present work Ey is adopted as the fundamental variable to be used to determine the 
QNM spectrum of polar electromagnetic perturbations. Having made this choice, one then 
writes equations in terms of the independent variable u = r^/r. With this, Eqs. ( p. 10 ) and 
(3.11) written for Ey lead to 



where Id and q are respectively the normalized frequency and wavenumber, defined by Eqs. 
(|3^), and \) is given by Eq. ( |l9|) . 

According to the Niihez-Starinets QNM definition [58], once one has found the funda- 
mental perturbation equations for a field on the AdS spacetime, the next step is establishing 
explicit relations between the perturbation variables and the corresponding retarded Green 
functions in the holographic GET. It is exactly from these relations that will emerge the 
boundary conditions to be imposed on the perturbation fields at infinity, viz, the conditions 
that lead to the singularities of the two-point correlation functions in the boundary field 
theory, and consequently to the quasinormal frequencies of the fluctuation modes. Such a 
task is performed in what follows. 

3.2 i2-current correlation functions 

In the case of electromagnetic perturbations, the AdS/GET correspondence [1-3] tells 
that, in the strong coupling, large limit, the information on the thermal correlation 
functions of the i?-symmetry currents are encoded into the electric field components Ej (j = 



x,y), which are solutions to the differential equations (3/7) and ( 3.14 ), respectively. It can 
be shown that, close to the horizon (n ~ 1), such functions are given approximately by Ej = 
j^±im/3^ where the negative (positive) exponent corresponds to ingoing (outgoing) waves. 
Moreover, depending on the sign of the exponent, the boundary values of the perturbation 
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functions act as sources of retarded or advanced Green functions in the dual CFT. To 
compute the retarded two-point functions, one has to opt for the negative exponent. It 
is also necessary to know the asymptotic form of the perturbation functions close to the 
infinite boundary (u ~ 0). A simple analysis shows that the solutions of equations ( |3.7| ) 
and (3.14) which satisfy incoming-wave condition at horizon present the following behavior 



around u = 0: 

= A^)(tt),q) + ... +%)(tt),q)M + (3.15) 
Ey = A(y){tti,q) + ... +B(^y){n),q)u + (3.16) 

where ellipses denote higher powers of u for each one of the independent solutions. Symbols 
^(j)(tD,q) and ;B(j)(tt), q), introduced in the above equations, stand for the connection 
coefficients associated to the corresponding differential equations for and Ey. 

To proceed further and calculate the correlation functions, the electromagnetic action 
at the boundary needs to be determined. It is usual to split the action as Sem = S^orizon + 
S boundary Uslug the equatious of motion and the preceding definitions it follows 



X ,. f dtvdq 
— iim 



^boundary ^ ^^^^^ J ^^^^ 

where 



E Ju, k)Ey{u, -k) H 2^x{'^^ k)Ej;{u, -k) 



SttTL^ (2iV)3/2r 



(3.17) 



is the electric susceptibility of the dual system [63,67]. 

In order to apply the Lorentzian AdS/CFT prescription of Ref. [37], the asymptotic 
solutions ( |3l5| ) and (CT ) are used to write the derivatives of the electric field in terms of 
the boundary values of the three-vector potential A^{k) = A^(u — > 0,A;). The i?-current 
correlation functions C^j/ are proportional to the coefficients of the terms containing the 
product Al{k)Al{-k) that appears in the action ( |3.17] ). It is then found: 



^ ^ q^ gfa)(tt),q) ^ tt)2 %)(w,q) 

" ^(tt)2-q2)^(^)(tt,,q)' yy ^(tt,2-q2)^(^)(rt,,q) 

„ __ tug gfa)(tT3,q) _ g(x)(i3,q) 

"-'y- ^(tr2-q2)^(^)(tt,,q)' "^^^ " ^^(,)(tt,, q) " 



(3.19) 



Moreover, for a (2 -|- l)-dimensional CFT at finite temperature, the current-current cor- 
relation functions can be written in terms of the transverse and longitudinal self-energies 
n-^(tD,q) and n^(rD,q), respectively (See Appendix ^ for a summary of such relations). 
Hence, comparing Eqs. ( p^ to Eqs. (^)-(Q of Appendix one finds 

ttT/^ a ^W(tt',q) . ^(y)(tt',q) 

^(^)(tt),qj ^(y)tD,q) 

These results show that the retarded two-point correlation functions are fully determined by 



- 12 - 



the ratio between the connection coefficients of equations (3/7) and ( 3.14| ). Furthermore, the 



poles of the thermal correlation functions are given by the zeros of the coefficients ^(a;)(tt), q) 
and ^(y)(tt),q). According to Ref. [58], the poles of C^jy define the electromagnetic QNM 
frequencies of the black hole localized in the AdS spacetime. Such frequencies are then 
obtained by imposing Dirichlet boundary conditions on the electric field components 
and Ey at u = 0, with E^ and Ey being functions that satisfy also an incoming-wave 
condition at the horizon. 

3.3 QNM and the gauge- invariant variables 

At this stage one could ask whether imposing Dirichlet conditions at the boundary 
(n = 0) on the Regge- Wheeler- Zerilli (RWZ) variables would produce the same 

QNM spectra as the spectra obtained by imposing the same boundary conditions onto 
the Kovtun-Starinets (KS) variables Ex^y For the transverse electromagnetic sector, the 
answer to this question is quite easy to find. In fact, variables and Ex are proportional 
to each other, so that both of the obtained QNM spectra, either using the KS variable or 
using the RWZ variable, are identical. In the case of polar electromagnetic perturbations, 
equations for the KS variable Ey and for the RWZ variable do not have the same 
form and, in addition, Ey and are independent variables, so that the answer is not 



immediate. In fact, as it is shown below (see Sects. 3.4 and p.5|), the QNM spectrum 



obtained from Ey is different from the QNM spectrum obtained from 

3.4 Dispersion relations for the hydrodynamic QNM 

The hydrodynamic limit of perturbations corresponds to the small frequency (tt) ^ 1) 
and small wavenumber (q ^ 1) region of the spectrum of the respective Fourier modes. In 
general, the quasinormal modes can be classified according the behavior of the dispersion 
relations in the hydrodynamic limit, and in this respect there are two classes. There is a 
set of QNM for which the frequency lt)(q) vanishes when q — > 0. Such modes are named 
here hydrodynamic quasinormal modes. But there is another kind of QNM for which the 
corresponding frequency in the long-wavelength limit is nonzero. To distinguish these two 
kind of modes from each other, the modes belonging to the later kind are denominated 
non-hydrodynamic quasinormal modes. In this section, the dispersion relations of the 
electromagnetic hydrodynamic QNM are studied by means of analytical and numerical 
methods. The electromagnetic non-hydrodynamic QNM shall be object of study in the 
next section. 

From the CFT point of view, it is expected that at least one of the electromagnetic 
QNM should show the typical behavior of a diffusion mode in the hydrodynamic limit. 
Trying to find such a mode, one then looks for solutions to Eqs. (3.7) and ( p. 14 ) in the 



form of power series in tr and q, under the assumption ro ~ q. Written in terms of the 
variables Fj = (for j = x,y), which are more appropriate for the present analysis, 

perturbation equations (|3.7| ) and ( |3.14D may be cast as 

2 

F- + — {2itv - Mj) Fj + ^ [itv{2u + - Saju^) + ^^{1- u^) - q^t)] Fj = 0, (3.21) 
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where = 1 and ay = lt)^/(rD^ — q^f)). After relabelling parameters as W ^ Xto and 
q ^ Aq with A <C 1, it is assumed that solutions of Eqs. ( |3.21| ) can be expanded in the 
form 

Fj{u) = F^{u) + XF^iu) + X'^Ffiu) + (3.22) 

where the coefficients F^{u), with a = 0, 1, 2, represent arbitrary functions of variable 
u, and which are also homogeneous functions of degree a on ttJ and q. 

The boundary condition of being ingoing waves at the horizon imposed on Ej, when 
translated to the new functions Fj, implies their dominant terms in expansion ( ^.22 ) must 



assume constant values close the horizon (u ~ 1). Then, in terms of the expansion ( 3.22| ) 
one has the following conditions 

F°(l) = constant, = if (1) = ... = 0. (3.23) 

It is now possible to solve Eqs. ( |3.21| ) order by order and, after imposing the boundary 
conditions given in Eqs. ( |3.23| ), the following expansions are found: 



Ey = ar""/' 



\/3 /tt l + 2u\ ixo ^ l + u + u^ 2s 
I ( arctan -I m ' 'T"-^' 



ivo — arctan — ^ H In h 0{xv^) 

3\3 y3/2 3 



(3.24) 



iq^ ,^ , . \/3 / vr 1 + 2u\ ivo ^ l + u + v? n 

1 + — l-u -?tt)— arctan — + — In + 0(^2 

W ^ ^ 3 V3 V3 y 2 3 ^ 

(3.25) 

where Cx and Cy are arbitrary normalization constants. One finds from Eq. ( |3.24| ) no 
solution satisfying the Dirichlet condition at the AdS spacetime boundary, i?x(0) = 0, and 
which is at the same time compatible with the hydrodynamic approximation tr, q ^ 1. This 
means there is no axial electromagnetic hydrodynamic QNM, and no i?-charge diffusion 
in the transverse direction to the spatial wave vector, as expected from the GET point of 
view. On the other hand, the condition Ey{0) = and Eq. (|3.25| ) lead to^ 

ro = -iq^ =^ uj = -^q\ (3.26) 

from where one can read the diffusion coefficient D = 2>/AttT. It is worth noticing that this 
diffusion mode is not found if one uses the RWZ master variable ^^^^ instead of Ey. 

As seen above, the hydrodynamic limit of perturbation equations comprises a very 
special interval in the space of parameters tt) and q. Besides the physical relevance of this 
regime, it corresponds to the very rare situations where analytical expressions can be found 
for the quasinormal frequencies. In the great majority of cases, numerical methods have 
to be employed in order to find the complete dispersion relations tt) x q. In the sequence, 
the Horowitz-Hubeny method [31] is used to compute the dispersion relation tr = —ivoi{q) 
for the electromagnetic hydrodynamic QNM, which in the limit of small wavenumbers 



^This result was also found through direct calculation of the hydrodynamic limit of correlation fun ction s 



by Herzog [63,64]. The comparison to that result is in fact a test for the analysis performed in Sect. 3.2 
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3q/47rT 

Figure 1: The dispersion relation for tlie only electromagnetic hydrodynamic QNM (solid line), 
which is purely damped, in = —ivoj, and corresponds to a polar perturbation. The dotted line 
is the diffusion mode ir/ = q^, which approaches the quasinormal frequency in the hydrodynamic 
limit ir, q ^ 1. 

corresponds to the diffusion mode found above. The result is shown in Fig. |l|. One sees 
the deviation of the exact dispersion relation curve (solid line) from the hydrodynamic limit 
curve tvj = (dotted line). Another interesting fact is that the quasinormal frequency 
tv{q) disappears for q larger than approximately 0.557. This is characteristic to all the 
electromagnetic purely damped modes, as analyzed in the next section (see Table [l|). 

3.5 Dispersion relations for the non- hydro dynamic QNM 

As mentioned earlier, one of the goals of the present analysis is to obtain the elec- 
tromagnetic quasinormal modes of plane-symmetric AdS4 spacetimes and to compare the 
present results with the results of Ref. [53]. As verified in the hydrodynamic limit (Sect. 
|3.4| ), the QNM spectra calculated here may be quite different from the spectra obtained 
in that work because of the use of different fundamental variables, and therefore a more 
careful search study on the dispersion relations of these modes is justified. 

3.5.1 Purely damped modes 

For small values of q, electromagnetic perturbations of AdS black holes present a 
special set of modes which are purely damped. These are not usual QNM since the real 
part of the frequencies vanishes eliminating the oscillatory behavior of the perturbations 
which is characteristic of QNM. Furthermore, the frequencies tt3(q) of such modes cannot 
be, in general, associated to hydrodynamic poles since most of the purely damped modes 
have nonvanishing frequencies in limit as the wavenumber q goes to zero. To distinguish 
from the regular QNM, the purely imaginary frequencies of both the axial and the polar 
perturbations shall be labelled by a special quantum number, n^, that assumes just integer 
values, starting by = for the hydrodynamic diffusion mode. As it was shown above, 
the axial sector of electromagnetic perturbations does not present quasinormal modes in 
the hydrodynamic limit, so that the set of axial purely damped modes starts at Ug = 1. 
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0.2 0.3 
3q/47tT 

Figure 2: The dispersion relations for polar (solid lines) and axial (dashed lines) purely damped 
electromagnetic QNM. The dotted (lowest) line is the diffusion mode Wg — (\^, which approaches the 
Us = quasinormal frequency in the hydrodynamic limit ir, q ^ 1. The insert shows the behavior 
of the dispersion relations for higher quasinormal frequencies. 



An interesting property of electromagnetic QNM of AdS black holes in (3 + l)-di- 
mensional spacetimes has been recently discovered by Herzog and collaborators [67]: The 
current-current correlators are analytical functions at q = 0, meaning that there are no 
quasinormal frequencies for null wavenumber. It can be shown that such a property is 
a consequence of the well known duality relation between electric and magnetic fields in 
vacuum. In fact, using the invariance of Maxwell equations under the duality operation, 
electric field <-> magnetic field, and the invariance of the correlation functions under rota- 
tions in the case of null wavenumber (zero momentum), it was shown that the transverse 
and longitudinal self-energies, n^(m,0) and n^(w,0), are well behaved functions of the 
frequency for all values of to. 

On the other hand, as verified through the numerical results for purely damped modes, 
there are quasinormal frequencies even for wavenumbers very close to zero. In fact, as shown 
in Fig. ^, the small wavenumber limit (q = e, with e very small but non-zero^) of the cor- 
responding purely imaginary quasinormal frequencies, tr = —itVs, is given approximately 

by 

3 

tfs = -^ris =^ ids = = '^ttTus, (3.27) 

where a;„^ are the Matsubara frequencies of a generic quantum bosonic system. Moreover, 
as it is also seen from Fig. ||, the hydrodynamic pole = is similar to other purely 
damped modes. As a matter of fact, the only property that distinguishes a particular 
purely damped mode from another is the behavior of these modes around q = 0: The only 
QNM satisfying the condition limq^ott'(l) = is the hydrodynamic mode. 

At the opposite side of the quasinormal spectrum, i.e., for larger values of q, there are 
saturation points at a maximum wavenumber value, qunn beyond which the specific mode 
disappears (See, however. Figs. ^ and|5|). This seems to happen for everyone of the modes, 

^In this specific case, the time of computation spent by the numerical code based on the Horowitz-Hubeny 
method [31] written to find the quasinormal frequencies becomes very large as e approaches zero. 
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Hum ^ J-u 


Itit ( iTifpi'val ) 
* ^ iiin \ ^^^^^^ * C'-L / 




qiim X 10^ 


roiim (interval) 


(0,1) 


557.319 


[0.648111, 0.648429] 


(1,2) 


339.330 


[2.04771, 2.04811] 


(2,3) 


162.034 


[3.52507, 3.52562] 


(3,4) 


71.8726 


[5.01701, 5.01788] 


(4,5) 


31.0102 


[6.51286, 6.51384] 


(5,6) 


13.1892 


[8.00994, 8.01169] 


(6,7) 


5.55917 


[9.50785, 9.51033] 


(7,8) 


2.32839 


[11.0057, 11.0100] 


(8,9) 


0.970660 


[12.5041, 12.5097] 


(9,10) 


0.403180 


[14.0009, 14.0114] 



Table 1: Approximate limiting values of frequencies ro = — ilfiim arid wavenumbers quj^ for purely 
damped electromagnetic modes. The brackets indicate that the actual limiting values lie between 
the two indicated endpoints in each case. 



with qiijn decreasing for higher overtones (cf. Fig. The curves representing dispersion 
relations associated to two different but contiguous modes meet each other exactly at the 
saturation point, i.e., the two dispersion relation curves coincide at that point. The axial 
modes group in pairs according to the relation = {(1, 2), (3, 4), (5, 6), ...}, while the polar 
modes are paired as Ug = {(0, 1), (2, 3), (4, 5), ...}. The approximate hmiting values toiim 
and qiim for = 0, 1, 2, 10 are shown in Table |l[ 

The existence of a meeting point between two contiguous dispersion relation curves 
suggests that for the special wavenumber values q = qum the corresponding quasinormal 
frequencies Id = — itriim represent double poles of the CFT current-current correlation 
functions. A strong support to such a conclusion comes from the behavior of the connection 
coefficients A(^j^{n,q) as a function of tr = — imj for small values of q and Idr = 0. 
In Fig. ^ it is shown the profile of A(^x)i^: l) for 1 = qum — 0.0132 and < ro/ < 9. For 
this wavenumber, there are six purely imaginary quasinormal frequencies encompassing the 
= 1 to ris = 6 axial QNM, which by definition are the points where ^(2,)(tt3,q) = 0. 
In particular, for = 5 and Ug = 6, the zeros of A(^x){^^^) coincide, indicating that the 
corresponding quasinormal frequencies are identical (see also Fig. |2|). 

The multiplicity of specific quasinormal frequencies as poles of the correlation functions 
may also be verified through the derivatives of ^^^(tt), q) with respect to to. This was done 
numerically, by taking fixed values of q and letting q qiim. It was found that the first 
derivative of ^(j)(tt),q) in relation to ro vanishes when tr = — itoiim and q = qiim, but 
it is not zero for other values of the wavenumber. This result proves that to = — itDum 
corresponds to, at least, a double zero of ^(j)(lt),q), and consequently, to a double pole of 
the corresponding current-current correlation functions. 

3.5.2 Ordinary quasinormal modes 

The electromagnetic perturbations of AdS black holes present also a family of regular 
(ordinary) quasinormal modes whose frequencies have nonzero real and imaginary parts. 
The numerical results for the quasinormal frequencies of the first five regular modes are 
shown respectively in Figs. |^ and ^ for axial and polar fluctuations. 

The form of the dispersion relations tt3/j(q) and tD/(q) indicates a connection between 
the electromagnetic ordinary QNM and the family of purely damped modes discussed in the 
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Figure 3: The connection coefficient ^(j.)(tT),q) for roj^ — 0, Xvi — (0,9) and q = q^j^ ~ 0.0132. 
The points represent the Us = 1, 2, .., 6 purely imaginary quasinornial frequencies associated to the 
axial electromagnetic perturbations. Note the coincidence of the points corresponding to rtg = 5 
and Us — &, indicating a possible doubleness of the related quasinormal frequency. 



last section. As shown in Figs. |^ and ^, each regular quasinormal frequency only appears 
for q larger than a minimum wavenumber value, which (to a good approximation) coincides 
with the limiting value of the corresponding pair of purely damped modes, qum. That is to 
say, all the dispersion relations, rD(q), for the ordinary QNM begin at the points (Wiim, qum)) 
with the real parts starting from zero value, 'cor{(\ Inm) = 0, while the imaginary parts 
start at tt3/(q lum) ~ '^Um- For higher wavenumber values, the ordinary electromagnetic 
modes show a sequence of quasinormal frequencies whose imaginary parts grow with the 
principal quantum number n. In this respect. Figs. ^ and |5| show a simiral behavior for 
the real parts of the quasinormal electromagnetic frequencies. 




0.5 1 1.5 2 0.5 1 1.5 2 



3q/4jtT 3q/4jtT 

Figure 4: The real (left) and imaginary (right) parts of the frequencies to = 3a'/47rT for the first 
five ordinary axial electromagnetic modes as a function of the normalized wavenumber q = iq/A-KT. 
The quantum number n arranges the regular polar QNM in increasing order of values of tn/. In the 
right it is also shown the frequencies VOs — SlOs/AttT associated to the axial electromagnetic purely 
damped QNM. 
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0.5 1 1.5 2 0.5 1 1.5 2 

3q/47[T 3q/47lT 



Figure 5: The real (left) and imaginary (right) parts of the frequencies to = Slu/AttT for the first 
five ordinary polar electromagnetic modes as a function of the normalized wavenumber q = 3g/47rT. 
As for the axial modes, the quantum number n arranges the regular QNM from lower to higher 
values of roj. In the right it is also shown the frequencies tOg = 3Us/4:-jtT associated to the polar 
electromagnetic purely damped QNM. 
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Figure 6: Graphs of the dispersion relations of the first five regular polar electromagnetic QNM 
for large values of q. The figure on the left hand side shows the real part of the frequency, tUfl, 
while the graph on the right hand side is for the imaginary part, tt)/. 



The quasinormal frequencies found here show that AdS black holes are not good oscil- 
lators. As it is well known, an interesting way of measuring the quality of an oscillator is by 
means of its quality factor Q = tt)ij/2tx)/. In general, in the region of small wavenumbers, 
the electromagnetic QNM have very small quality factors, Q <C 1. For instance, taking 
q = 0.557319 and considering the fundamental polar mode one obtains Q = 7.24 x 10~^, 
a quality factor typical to highly damped oscillatory systems. On the other hand, quality 
factors of the order of unity are found for large wavenumber values, such as Q = 1 for q 
around 1.11, and Q = 85.9 for q = 40. 

Prom the holographic field theory point of view, the real part of the frequencies may be 
interpreted as quasiparticle excitation energies in the dual plasma defined at the conformal 
boundary of the AdS spacetime. However, such an interpretation only makes sense for 
excitations, or quasinormal modes, with large quality factors. In fact, according to Hcisen- 
berg uncertainty principle, the uncertainty in the energy of a quasiparticle is of the order 
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of uii (in units of h). Hence, quality factors smaller than unity imply in energy uncertain- 
ties larger than the energies of the quasiparticles themselves, making the interpretation of 
energy excitations as quasiparticles meaningless. 

In the Q ^ 1 regime, where the quasiparticle interpretation is feasible. Fig. ^ shows 
that the dispersion relations of the ordinary electromagnetic QNM frequencies have real 
parts approaching straight lines of the form vor = (\ + ^n, where b„ depends only on the 
specific mode n. This means the functions ujii{q) approach the usual energy-momentum 
relation associated to a zero rest mass particle, u^r = as T ^ 0. Furthermore, the 
characteristic damping time (r = of the electromagnetic fluctuations diverges in the 

limit q — > oo, i.e., the functions ujj{q) tend to zero for large wavenumbers. All of the 
above results are consistent with the expected properties of poles of correlation functions 
in quantum field theories at zero temperature [67]. 

4. Gravitational quasinormal modes 

Even though gravitational perturbations of plane-symmetric black holes in (3 -|- 1)- 
dimensional AdS spacetimes have been analyzed in some extent [53-55,63,64], the quasi- 
normal-mode dispersion relations, with the use of the KS variables, were not found yet, 
and hence the comparison with the spectra obtained by using the RWZ gauge-invariant 
variables was not performed. This is done next. 

4.1 Fundamental equations for gravitational fluctuations 

As usual, gravitational perturbations are described here in terms of linear metric fluctu- 
ations, which means the metric for the perturbed spacetime is written as gMN = guN^^MN, 
where Hmn is considered as a perturbation in the background metric g"fjv) given by Eq. 
( |2.6| ). Among the variety of possible gauge choices in studying metric perturbations, an 
interesting choice is the so-called radial gauge, in which the coordinate system is chosen in 
such a way that hrt = hrx = hyy = h^r = 0. Since one of the aims here is to investigate the 
relation among the perturbations of plane-symmetric AdS4 black holes and the correlation 
functions in the dual CFT, it is convenient to use the radial gauge formalism to study the 
metric fluctuations. A brief description of this formalism is given in the following. 

As it was done in the previous section when studying the electromagnetic perturba- 
tions, the isometrics of spacetime (|2.6| ) allow Fourier transforming coordinates t, x and y, 
and writing metric fluctuations /imjv as 



Again the wave three- vector may be chosen as = (— w,0,g), and hence metric per- 
turbations can be split into two disjoint sets. Namely, the axial (transverse) sector of 
gravitational perturbations is characterized by the quantities htx and hy^, and the polar 
(longitudinal) sector of perturbations is composed by htt, hxx, hyy, and hty 





(4.1) 
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4.1.1 Axial perturbations in the radial gauge 

Now one needs to find evolution equations for the axial (transverse) gravitational per- 
turbations in the radial gauge, that are composed by the metric fluctuations htx and hyx- 
The task is carried out following the usual procedure of the theory of linear perturbations. 
The linearized Einstein equations corresponding to the axial sector of gravitational pertur- 
bations yield a set of three coupled differential equations for htx and hyx [63] . Of course, 
such equations are not independent from each other, and one of them may be eliminated as 
a combination of the other two. A resulting system of linearly independent equations (after 
Fourier transforming them) which is interesting for the present study is the following: 

Htx + ^Hy, = 0, (4.2) 
< - - ^ {tvHy, + qHtx) = 0, (4.3) 

where Htx and Hy^ are defined by 

Htx htx^ Hyx hyx- (4'4) 

As in the case of electromagnetic perturbations, a mandatory condition that a candi- 
date for fundamental variable must satisfy is being gauge invariant, which for metric per- 
turbations means the candidate has to be invariant under infinitesimal coordinate transfor- 
mations. Inspired once again in the work by Kovtun and Starinets [60], among the different 
combinations of axial functions Htx and Hyx which furnish gauge-invariant quantities, one 
takes 

Zi = i {qHtx + ojHyx) (4.5) 

as the fundamental gauge-invariant function of the axial gravitational perturbations. 

Decoupling the system of differential equations (^]^) and (l4.3| ) in terms of the funda- 
mental variable Zi, it results the solely second-order differential equation 

nD(m^-q^D) + = ^'-'^ 

Solutions to this equation satisfying the QNM boundary conditions are studied below. 

4.1.2 Polar perturbations in the radial gauge 

The polar (longitudinal) sector of gravitational perturbations in the radial gauge is 
described by metric fluctuations htt, hxx, hyy and hty. These components of the metric 
perturbation tensor are used to deflne new quantities 



1 r2 r2 7-2 

-h TJ -—h TJ -—h TJ - — 

J. 'Hti -'^xx — ^2 "'xxi ^yy — ^2 TO' ty — ^2 



which are more appropriate to deal with during calculations to obtain perturbation equa- 
tions. Hence, the polar components of linearized Einstein equations furnish a set of seven 
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coupled equations for the variables defined in Eqs. (4/?). Only four of such a set are in- 
dependent equations and, among the possible choices, the more interesting for the present 
work are 



2nn3q , ua(u) , , 4rDfi 

h[u) 2qf)6(ii) ^ qK^) 



[4ffi2i?,, + c(n)//it] - (tt)i?,, + XoHyy + 2qHty) , (4.10) 



(4.11) 

where, as above, the primes denote derivatives with respect to the variable u = r^/r, and 
coefficients a{u) = 3u^ — 12u — Ato"^, h{u) = + 3, and c{u) = Ato'^ — (\^b{u) were introduced 
to simplify notation. 

Finally, a gauge-invariant function Z2 is built as a particular combination of the metric 
perturbations 

Z2 = AuoqHty + 2uJ^Hyy + [q\3 -I))- 2uj^] H,, + 2q^l)Hu, (4.12) 

for which, uncoupling the equations of motion ( [4.8D -( pyTl| ), it is found the following second- 
order differential equation 

„ 4tt)2(2 + u3) + q2d(n) ^, , 4tt)" + q4[)fe(n) - q^n) ^ „ 
ux)c(u) r)"^c(u) 

where d{u) = Av? — bu^ — 8 and e{u) = 9n'^f) -|- tt)2(8 — 5n^). Equations ( ^^ and ( |4.13D are 
the fundamental equations that are going to be used in the next sections to compute the 
QNM spectra associated to gravitational perturbations of plane AdS4 black holes. 



4.2 Stress-energy tensor correlation functions 

For the gravitational perturbations, the AdS/CFT correspondence establishes a re- 
lation among the solutions of Eqs. ( [4. 61) and (|4.13D and the stress-energy tensor of the 
dual GET. From this relation, the stress-energy tensor correlators can be determined, and 
in order to do that the explicit form of the fields in the bulk AdS spacetime has to be 
known. More precisely, in order to impose the ingoing-wave condition at the horizon, and 
to map AdS to GET quantities at the boundary of the spacetime, the asymptotic form of 
the metric perturbation functions close to the horizon and at the boundary are necessary. 

In the horizon neighborhood (tt ~ 1), the gravitational gauge-invariant variables Zi 
and Z2 have a similar behavior as the electric field components (see Sect. ^^), viz, Zi^2 ~ 
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j^±iTO/3 ^j^g electromagnetic case, to compute the retarded Green functions, one has 

to choose the solutions corresponding to the negative imaginary power, Zi^2 ~ 
the other hand, at the conformal boundary of the AdS spacetime, the metric fluctuations 
are such that 



Zi =A(^i){tv,q) + ... +i3(i)(tt),q)u3 + 
Z2 = A(2){rv, q) + ... + %)(ID, q)u^ + ... , 



(4.14) 
(4.15) 



where ellipses denote higher powers of u, and quantities ^(i)(rD,q) and ;S(i)(tt3, q), and 
-^(2)(rt') q) and B(^2)i^^ 1) a-re the connection coefficients related to the differential equations 
( |4.6| ) and ([4.13 ), respectively. 

For the remaining of this section, as usual, the gravitational perturbations are split 
into axial and polar sectors and the analysis of the corresponding actions, coming from Eq. 
( |2.9| ), are performed separately for both of the perturbation types. 



4.2.1 Axial sector 

It is well known that in the calculation of two-point correlation functions from the 
gravitational action only quadratic terms in metric perturbations need to be considered. 
Moreover, according to the Lorentzian AdS/CFT prescription [37] (see also Ref. [80]), in 
order to obtain the CFT retarded correlators the relevant terms are the quadratic terms 
in the derivatives of Hnu- Hence, collecting all of the contributions coming from the grav- 



itational part of action (2.9), one gets 



5(2) 



du d^X—p; 



tx 



+ .... 



(4.16) 



where 



P 



81 



(4.17) 



is interpreted as the pressure of the dual plasma [64] . 

Now expressing functions and Hy^ in terms of the axial fundamental variable Zi 
through Eqs. (|4.2| )- (|4.5| ), substituting the resulting relations into Eq. ( [4.16| ), and making 
use of the fundamental equation ([4.6| ) , it is found the following action (at the boundary) 



0(2) 

'^boundary 



— lim 

2 n^O 



dto dq 



Zi{u, k)Zi{u, —k) + 5, 



(2) 



(4.18) 



where the contact terms represented by S^t do not carry derivatives of the metric perturba- 
tion functions. In the calculation of the correlation functions, after Fourier transformation, 
the contact terms give rise to derivatives of the Dirac delta function. Their removal can be 
done through the holographic renormalization, with the inclusion of appropriate counter 
terms in the supergravity action [81] . 
Besides using Eq. 



), the asymptotic expansion given by Eq. ( 4.14 ) is used to write 
the derivative of the gauge-invariant quantity Zi in terms of the boundary values of the 
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perturbation fields H^j^{k) = Hf^i^{u — > 0, fc), and then the AdS/CFT prescription [37] can 
be apphed to the present case in order to calculate the retarded correlation functions of 
the holographic stress-energy tensor T^'^. The result is 



a 



tx,tx 



-3P 



^(i)(tt',q) 



G. 



tx,yx 



3P 



G. 



yx,yx 



(tr2-q2)^(^)(tt,,q)' 

tug g(i)(tD,q) 
(lt)2-q2)^(^)(tt,,q)' 

tD^ g(i)(tt),q) 
(tt)2 -q2)^(i)(lt,,q)' 



-3P 



(4.19) 
(4.20) 
(4.21) 



As it happens for the current-current correlations functions, one can find general ex- 
pressions for the two-point thermal functions associated to the stress-energy tensor which 
hold for any scale invariant (2 -|- l)-dimensional field theory (see Appendix^. For fluctu- 
ations of the transverse momentum density in the CFT one has the correlators 



Gtx,tx 

G. 



2 tt)^ 



roq 



^tx,yx 



2(tr2 -q2) 



■G'i(lD,q), 
Gi(ti3,q), 



yx,yx 



2(tt)2-q2; 



Gi(lD,q). 



(4.22) 
(4.23) 

(4.24) 



Therefore, by comparing the general expressions of Eqs. ( [4.22| )-(4.24) to the results given 
in Eqs. ( 4.19D -( 4.21| ), the following scalar function is found 



Gi(tD,q) = -6P 



^(i)(tt),q) 



(4.25) 



It is then seen that Dirichlet condition imposed on the fundamental variable Zi at the 
boundary, Zi{0) = A(^i^{to,i:\) = 0, leads straightforwardly to the poles of the correlation 
functions Gtx,tx, Gtx,yx and Gyx,yx- As a consequence of this result, such a requirement also 
yields the quasinormal spectrum associated to the axial gravitational perturbation modes 
of plane-symmetric black holes. 



4.2.2 Polar sector 

The procedure to be applied to the polar sector of gravitational perturbations is the 
same as for the axial sector. The starting point here is the part of the boundary gravi- 
tational action built with the quadratic terms in the polar metric perturbations, which is 
given by [64] 



0(2) 

^boundary 



P 

— lim / (Px 

2 u^o 



^ — '^Hiy + HttHxx + HuH. 



yyi 



^ {Hxx Hyy) 



{^ly + ^XxHyy — HtfHxx " HftHyy^ 



(4.26) 
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By using the relation among polar metric fluctuations and the gauge- invariant variable Z2, 
Eq. ( [4.12 ), and the equations of motion (4.8)-(4.11), the boundary action ( 4.26| ) can be 



cast into the form 

Siounclary = ^ 1™ / ^— — 3 Z'^iu, k)Z2{u, - k) + S^cl , (4-27) 

where the contact terms S^l do not contain derivatives of metric perturbations. One now 
uses the asymptotic expansion ( |4.15| ) to write the derivative of the gauge-invariant variable 
Z2 in terms of boundary values of the polar metric perturbations H^^ik) . After substituting 



the resulting expression into the action ( 4.27 ), the appropriate functional derivatives^ of 
the action with respect to the independent fields H^^{k), H^y^k), Hyy{k) and H^^{k) are 
performed. Notice, however, that in the case of polar gravitational perturbations, the 
use of the AdS/CFT prescription is not direct. In fact, it is first necessary to identify 
explicitly how the metric perturbations couple to the stress-energy tensor at the boundary. 
As discussed in Refs. [82,83], such a coupling is given by 

-\jdt<fx h'^^T'^, = -ljdtcPx [HlT'' + + H^yTyy + 2H%T'y] . (4.28) 

Taking this coupling into account, the covariant components of the polar correlation func- 
tions are found 

G^v^ap = Q^lv,a|3G2{vo,(:\), (4.29) 
where the scalar function G2(tt), q) is given by 

G2(lD, q) = -6p4^^4-^ + contact terms, (4.30) 

A2)("''l) 

and tensor Qf^u,a(5 is given in Appendix ^ (see also [60]). This result shows definitely that 
Dirichlet boundary condition imposed on the fundamental variable Z2 at infinity leads to 
the poles of the function G2(tt),q), and, by definition, to the quasinormal frequencies of 
polar gravitational vibration modes. 

4.3 QNM and the gauge- invariant variables 

As in the study of electromagnetic quasinormal modes, a comparison between results 
found using RWZ variables with the ones obtained using KS gauge-invariant variables 
to describe gravitational perturbations deserves to be made. As in that case, while the 
axial quasinormal spectrum is independent of the choice of the fundamental variable, polar 
quasinormal spectrum strongly depends on it. 

A strong evidence that the QNM spectra obtained using either RWZ or KS variables 
to describe axial gravitational perturbations are identical comes from the study of the 
hydrodynamic limit of such perturbations as performed in Ref. [54] and in Sect. (see 
below). As a matter of fact, even though different methods have been employed in each 



^Functional derivatives in the sense defined by the Lorentzian AdS/CFT prescription of Ref. [37]. 
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case, both of the quasinormal spectra present a typical hydro dynamic shear mode with 
diffusion coefficient D = 1/4ttT, independently if one uses variable Z'~\ or variable Zi. 
Furthermore, it can be shown that the explicit relation between RWZ and KS variables is 

Zi = ^dr [rZ^-^] , (4.31) 

and at the AdS spacetime boundary, axial fundamental variables and Zi are propor- 
tional to each other, 

^i(^)L=o = i2^"NL=o' (4-32) 

what proves that the two spectra obtained from and Zi are indeed identical to each 
other. 

The situation in the polar sector is quite diverse from what happens in the axial sector. 



As it is shown in the sequence of the present work (see Sects. and f4.5| ), RWZ and KS 



variables with the same boundary conditions generate different quasinormal frequencies. 



In particular, it is shown in Sect. 4.4 that the hydrodynamic limit of Z2 contains a sound 



wave mode which is not seen in the quasinormal spectrum obtained from Z'-^K 



4.4 Dispersion relations for the hydrodynamic QNM 



For the fluctuations of the stress-energy tensor in the dual field theory, hydrodynamics 
predicts a shear mode in the transverse (axial) sector and a sound wave mode in the 
longitudinal (polar) sector. As it is going to be shown below, these modes also appear 
in the gravitational QNM spectra of the plane black holes as long as one investigates the 
regime of small frequencies and wavenumbers, tD — > and q — > 0. Proceeding analogously 



to the case of electromagnetic perturbations studied in Sect. |3.4 the following change of 
variables is done 



u 



J = 1,2. 



(4.33) 



Functions Hi and H2 are then expanded in power series of Id and q. Besides being ap- 
proximate solutions to Eqs. (4^) and ([4.13| ) in the hydrodynamic limit, such series must 
represent ingoing waves at the horizon, namely, Hj(u)\^_ -, = constant. These conditions 
are fulfilled by the following expansions: 



\u=l 



Z2 = C2i) 



-ivo/3 



u 



4ilD(] 



+ 0(lt)2 



(4.34) 



(4.35) 



q^ 3 

where Ci and C2 are arbitrary normalization constants. Dirichlet conditions at the space- 
time boundary, u = 0, are now imposed on both the axial and the polar fundamental 
variables, Zi{0) = and Z2{0) = 0. The first variable leads to an axial quasinormal mode 
identified to the shear mode with 



UJ 



-iDq^ 



(4.36) 
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while the second variable furnishes a dispersion relation characteristic to a sound wave 
mode 

= ±-^g - ^ , (4.37) 
where D is the diffusion coefficient, given by 

D = = t4^, (4.38) 

e + P AttT' ^ ' 

with r\ and e being respectively the shear coefficient and the energy density of the dual 
system. 



hydrodynamic QNM 
shear mode 
algebraically special mode 



2 - 




Figure 7: The dispersion relation (imaginary part) for the axial gravitational hydrodynamic mode 
(solid line), in comparison to the shear mode, q^/3 (dotted line), and to the algebraically special 
mode, q^/6 (dashed Une). 



Among other interesting results that can be found, combining Eq. (|4.38| ) to the ther- 
modynamic Euler relation P = —e + Ts =^ s = (e + P) /T, one gets the ratio between the 
shear coefficient and the entropy density of the dual plasma, 



1 

s 



7] T 

e + P 



1 

4^' 



(4.39) 



or, in conventional units, r]/s = h/Ankg. This ratio is the same for all finite temperature 
field theories with a dual gravitational description in the AdS spacetime [84]. It is also 
speculated that rj/s = I/Att represents a lower bound -the KSS bound- of such a ratio for 
all fluids in nature. 

The complete dispersion relations for the gravitational hydrodynamic QNM are ob- 
tained by means of the Horowitz-Hub eny method [31] , and, for the sake of simplicity, the 
numerical results for the axial and polar sectors are analyzed separately. 



(i) Axial modes: The extension of the dispersion relation of the axial hydrodynamic QNM 
for large values of q was already done in Ref. [54], but for completeness it is also shown 
in Fig. ^. As the normalized wavenumber q reaches values beyond of the hydrodynamic 
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regime, the magnitude of the QNM frequency tv = —itvi increases faster than the mag- 
nitude of the shear mode frequency tv = — and for very large wavenumber values, 
q » 1, the hydrodynamic QNM frequency approaches the algebraically special frequency 
ro = —iq'^/6. 



10 



hydrodynamic QNM 
sound wave mode 




Figure 8: Real and imaginary parts, respectively, of the dispersion relation iTi(q) for the hydro- 
dynamic longitudinal gravitational mode (solid li nes). The dashed lines in both of the figures 
correspond to the sound wave mode, given by Eq. ( 4.37 ), and the dotted curve in the figure on the 
left hand side is the relation Wr = q. 
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Figure 9: The group velocity Cg = dtxiji/dq and the derivative dWi /dq as a function of the normal- 
ized wavenumber q for the polar hydrodynamic gravitational mode (solid lines). The dashed line 
in the figure on the left is the sound velocity for a (2-|-l)-dimensional CFT, Cg = 1/^2, and the 
dotted line in the same figure represents the speed of light, c = 1. 



(a) Polar modes: In this sector of perturbations, RWZ and KS variables submitted to 
the same boundary conditions generate different quasinormal frequencies. As seen above, 
the hydrodynamic limit of Z2 presents a sound wave mode which is not present in the 
quasinormal spectrum obtained from Z'-'^\ The extended dispersion relations for the lon- 
gitudinal hydrodynamic QNM are shown in Fig. The real part of the frequency clearly 
shows the transition from the hydrodynamic regime = q/V^, at low wavenumbers, to 
a regime characterized by collisionless dual plasma in which the dispersion relation tt)ij(q) 
approaches the ultra-relativistic relation vor = q. Between these two extreme regimes, the 
group velocity, defined by Cg = dtoji/dq, assumes values that are higher than the speed of 
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light. The graphs in Fig. ^ show that Cg > 1 for all wavenumbers larger than q ~ 1.336, 
and that the minimum decaying time (the maximum of tr/) corresponds to q ~ 3.213. 
Notice, however, that dtvu/dq surpasses the speed of light at wavenumber values lying 
outside the hydrodynamic regime and therefore that superluminal group velocity cannot 
be interpreted as the sound velocity in the corresponding media. 

4.5 Dispersion relations for the non- hydro dynamic QNM 

Conventionally, a non-hydrodynamic QNM is every mode for which the dispersion 
relation presents a gap in the limit q — > 0. That is to say, the quasinormal frequency lt)(q) 
tends to a nonzero value in the limit where the wavenumber goes to zero. These kind of 
gravitational perturbation modes are studied in this section. 

As done in the case of electromagnetic perturbations, the QNM obtained by using 
RWZ gravitational variables Z'-'^\ obeying Schrodinger-like equations, are compared to the 
QNM obtained by using the KS gauge-invariant quantities Zi^2, which lead to the poles of 
stress-energy tensor correlators in the M = 8 super- Yang-Mills field theory. In particular, 
the results obtained here from Zi^2 sue compared to the results of Refs. [53,54]. 

4.5.1 Purely damped modes 

The spectra of gravitational perturbations of plane-symmetric AdS/^ black holes do 
not present non-hydrodynamic quasinormal frequencies with vanishing real part. In fact, 
the only purely damped mode of the gravitational perturbations is the axial hydrodynamic 
QNM which was already investigated in the last section. 

4.5.2 Ordinary quasinormal modes 

As usual, the study of the dispersion relations of regular non-hydrodynamic gravita- 
tional QNM is more conveniently performed by considering axial and polar sectors of such 
perturbations separately. 

(i) Axial modes: As shown above, in the case of gravitational axial modes both of the 
fundamental variables Zi and Z^~^ yield the same QNM spectrum. Even though a detailed 
study of the fundamental non-hydrodynamic QNM (based on the RWZ master variable) 
was performed in Ref. [54], higher overtones of axial gravitational QNM were not fully 
investigated. Hence, the aim here is to complete the analysis by including such higher 
overtones. Fig. |l^ shows the numerical results for the dispersion relations of the first five 
axial quasinormal modes: n = 1,..., 5. The general forms of the curves are approximately 
the same for all values of n: At intermediate values of q, there is a local minimum in the real 
part of the frequency rDij(q), and for large values of the wavenumber, every dispersion rela- 
tion tT)R(q) tends to some straight line parallel to the ultra-relativistic energy-momentum 
relation, xor = q. 

(ii) Polar modes: Contrary to the axial perturbations, the longitudinal gravitational fluc- 
tuations present different spectra as one takes Z2 or as fundamental variable. The 
numerical results for these modes, based on the gauge-invariant variable Z2, are shown 



-29- 



14 



12 :_...?.= 5. 




' ' ' ' ' ' ' ' ' ' 1 I ' ' ' ' ' ' ' ' ' 1 

0123456789 10 0123456789 10 

3q/4jtT 3q/4jtT 



Figure 10: The first five quasinormal frequencies of non-hydrodynamic axial gravitational modes, 
ro = (3u;ij/47rr) — z(3cj//47rT), as a function of the normalized wavenumber q — 'iq/AnT. The 
quantum number n arranges the modes in growing order according to the strength of the imaginary 
parts of the frequencies. 
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Figure 11: The first five quasinormal frequencies of non-hydrodynamic polar gravitational modes, 
tn = {"iiOB./ ^T^T) — i(3w//47rT), as a function of the normalized wavenumber q = iq/AnT, ordered 
as in the case of axial gravitational modes. 



in Fig. |11|. Notice that now the real part of the quasinormal frequency Vor is a mono- 
tonic increasing function of q, showing neither local maxima nor local minima which, on 
contrary, appear in the quasinormal spectrum for the RWZ master variable Z'-^^ used in 
Refs. [53,54]. A comparison between the QNM spectrum found from the RWZ variable (cf. 
Ref. [54]) and that found from the KS variable used here can also be done through the data 
presented in Table |2[ It is clearly seen the similarity between the two spectra in the two 
asymptotic regions of wavenumber values. Both of the spectra are approximately the same 
for small values and also for large values of q. The main differences happen in the regime 
where the normalized wavenumber q is of the order of unity, which means that q is of the 
order of the blackhole Hawking temperature. Moreover, the real parts of the quasinormal 
frequencies in both of the spectra are relatively closer to each other when compared to 
the corresponding imaginary parts. This implies that, even though the QNM oscillation 
frequencies are essentially the same for both choices of variables, the decaying timescales 
(r = l/w/) are significantly smaller for the KS choice. 
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Kovtun-Starinets Regge-Wheeler-Zerilli 



q 




tvi 






0.004 


1.84942 


2.66385 


1.84945 


2.66384 


0.04 


1.84964 


2.66379 


1.85027 


2.66248 


0.4 


1.87207 


2.65770 


1.92488 


2.52658 


1 


2.00603 


2.62917 


2.03016 


1.92213 


2 


2.60256 


2.56803 


2.30526 


1.55218 


5 


5.57791 


2.25854 


5.25618 


1.27974 


10 


10.5703 


1.94304 


10.2839 


1.07342 



Table 2: Numerical results for the first non-hydrodynamic quasinormal mode associated to the 
polar gravitational perturbations. The second and third columns show respectively the values of 
real and imaginary parts of the frequencies obtained by using the KS variable Z2, and the last 
two columns present the results obtained by using the RWZ variable which are also shown in 
Ref. [54]. 

5. Final comments and conclusion 

One of the important issues dealt with in the present work is related to the choice 
of appropriate variables in order to determine the QNM spectra of AdS black holes. It 
is argued that Kovtun-Starinets gauge-invariant quantities, together with incoming- wave 
condition at horizon and Dirichlet boundary condition at infinity, should be used in order 
to find the correct quasinormal dispersion relations. The resulting spectrum for a given 
perturbation is in general different from what is obtained using other kind of variables with 
the same boundary conditions. 

In the case of electromagnetic perturbations, contrary to what is obtained using the 
so-called Regge-Wheeler-Zerilli quantities where the spectra for both the axial and 
the polar modes are the same, Kovtun-Starinets variables present a different spectrum for 
each mode type. Also, both these sectors of electromagnetic perturbations present purely 
damped modes whose dispersion relations, in the limit of small wavenumbers, approach the 
bosonic Matsubara frequencies co = —2iTiTns. This result can be compared to the "quasi- 
Matsubara" frequencies, a; = 27rTns(l — z), that have been found for zero wavenumber 
fluctuations of (4 -|- l)-dimensional black branes [58]. The studies show the emergence of 
infinite sequences of bosonic Matsubara frequencies for both (3-|-l)- and (4-|-l)-dimensional 
black branes, but the particular behavior of the QNM dispersion relations at zero wavenum- 
ber is quite different in each case. In (3 -|- 1) dimensions the real part of the frequencies 
is zero for very small wavenumber values, while in (4 + 1) dimensions real and imaginary 
parts of the frequencies are both finite at zero wavenumber. Moreover, as pointed out in 
Sect, ^is] and discussed in detail in Ref. [67], the invariance of Maxwell equations under the 
electric field <-> magnetic field duality operation in (3 -|- l)-dimensional spacetimes implies 
that there are no electromagnetic QNM at zero wavenumber. Such a duality invariance 
does not hold in higher dimensional spacetimes, what justifies the different behavior of zero 
wavenumber electromagnetic QNM found here when compared to the results of Ref. [58] . 

Other special property of electromagnetic fluctuations deserving to be mentioned here 
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is the cutoff in tlie dispersion relations of purely damped QNM at a particular value of 
the wavenumber, q = qum- This cutoff implies in an abrupt change in the behavior of the 
fundamental quasinormal mode, but not in the thermalization time r. For wavenumbers in 
the interval < q < qiim, parameter r is given by the first purely damped mode, r = 1/ws, 
while for wavenumber values above qiim; the characteristic decaying time are governed by 
the fundamental ordinary QNM, r = l/uoi. Since the imaginary parts of the frequencies 
of these modes are equal for q = qunn the thermalization time changes continuously for 
wavenumbers close to qum. 

The numerical results for the gravitational QNM show that the thermalization time of 
axial modes for wavenumbers in the interval < q < 1.935 is determined by the hydrody- 
namic mode. This means that, at least in the limit q ^ 1, where the shear mode is a good 
approximation for the QNM, the thermalization time r is a linear function of Hawking 
temperature, r ~ AirT/q^ . As in the case of electromagnetic fluctuations, the transition to 
the regime where the thermalization time is determined by the first regular axial QNM is 
continuous. Such a transition happens for the values q ~ 1.935 and ft)/ ~ 2.296. Also, at 
this point the thermalization time reaches its minimum value, r ~ 0.104/r. On the other 
hand, for polar gravitational perturbations the decaying time r is always determined by 
the hydrodynamic QNM which reduces to the sound wave mode in the small wavenumber 
limit q <C 1. 

Finally, the behavior of the group velocity Cs shown in Fig. ^, which is greater than 
unity for all q > 1.336, deserves further analysis. First note that apparent superluminal 
propagation of this type, which at first sight seems to violate causality, has been found in 
other relativistic quantum systems. For instance, Scharnhorst [85] has shown in the case 
of Casimir effect that when vacuum fluctuations obey periodic boundary conditions, the 
two-loop corrections to the polarization tensor lead to superluminal photon propagation. 
Also, it is argued that the physically meaningful propagation velocity and, consequently, 
the one that defines the light cones in spacetime, is the front wave speed w„f, which is given 
by the limit of the phase speed Uph = uj/q when a; ^ oo (see, e.g., Ref. [86] for a review). 
The graph on the left in Fig. ^ shows that Vph approaches c = 1 at high frequencies, and, 
therefore, the oscillations in the super- Yang-Mills plasma do not violate causality. 
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A. Correlators in a (2 + l)-dimensional CFT 

The Lorentz index structure of retarded Green functions of conserved currents and 
stress-energy tensor in a D-dimensional relativistic quantum field theory was discussed in 
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details in Rcf. [60] . It was shown that the field theory correlators can be expressed in terms 
of a set of scalar functions. For the present work, a particularly interesting example is that 
of a (2 + l)-dimensional finite temperature conformal field theory. In this specific case, 
with the wave three- vector in the form A;^ = (— a;,0, g), the transverse component of the 
current-current correlation functions can be written as 

C,i,i(fe) = n^(u;,g). (A.l) 

The longitudinal components in turn are given by 

Ctt{k) = j^^^^Il^{u,q), (A.2) 

CtAk) = --^^Tl\u,ql (A.3) 
(o;^ — q ) 

C^2Ak) = -^ ^U'^{u,q), (A.4) 

[to - q ) 

where Il^{uj,q) and n^(aj,g) are two independent scalar functions. All the correlators of 
transverse momentum density are expressed in terms of a scalar function Gi{u,q): 

1 g2 
z [LO q ) 

1 o;^ 

Gx^x^x^x^ik) = -^-i-^—^Gi{uj,q). (A.7) 

In a similar way, the correlators of longitudinal momentum density, energy density, and 
diagonal stress are determined by another scalar function G2{u>,q), 

Gnu,af3{k) = Qfiu,a(3{i^, 9)^2(0;, q), (A.8) 
where the components of the projector Qijtu,ap are given by: 

Qtt,tt = 2 (^2 _ g2)2 ' Qtt,tx^ = - 2 (a;2 _ g2)2 ' (^"9) 

n - n 

Qx^x^ ,x^x^ 2' Qx^x^,tx^ 2 ^^2 g2^ ' (-^•-'--'-) 



1 CJ^ ]_ ^4 

Qxia;i,x2a;2 = ~ 2 (w^ — (/^) ' Qx^x^,x'^x^ — 2 ((^2 _ (^2"j2 ' (-^•-'-^) 

1 a;^^' 1 w^g^ 



Qx2a;2,tx2 2 (^;2 _ ^2)2 ' Qtx^,tx2 ^ [J^ — g2-)2 ' (A. 13) 
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